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This paper discusses a new variable equation of state parameter leading to exact solutions of the 
Einstein field equations describing traversable wormholes. In addition to generalizing the notion of 
phantom energy, the equation of state generates a mathematical model that combines the generalized 
phantom energy and the generalized Chaplygin gas models. 
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I. INTRODUCTION 

Traversable wormholes, whose existence was first conjec- 
tured by Morris and Thorne in 1988 [![, may be defined 
as handles or tunnels in the spacetime topology linking 
different universes or widely separated regions of our own 
Universe. A renewed interest in wormholes is due in part 
to the discovery that our Universe is undergoing an accel- 
erated expansion (HQ, that is, a > in the Friedmann 
equation a/a = (— 47r/3)(p -I- 3p), using units in which 
c = G = 1. The cause of the acceleration is a negative 
pressure dark energy, a form of matter whose equation 
of state is p = —Kp, p > 0, and K a constant; p is the 
spatially homogeneous pressure and p the energy density. 
A value of K > 1/3 is required for an accelerated expan- 
sion. The case 1/3 < K < 1 is referred to as quintessence, 
while K = 1 corresponds to a cosmological constant |4(. 
Of particular interest is the case K > 1, referred to as 
phantom energy, since it leads to a violation of the null 
energy condition, an essential requirement for maintain- 
ing a wormhole p]]. Matter that violates the null energy 
condition is usually called exotic. Since the notion of dark 
or phantom energy ordinarily applies only to a homoge- 
neous distribution of matter, while wormhole spacetimes 
are necessarily inhomogeneous, phantom energy does not 
automatically qualify as a candidate for exotic matter. 
It is shown in Ref. [5], however, that an extension to 
spherically symmetric inhomogeneous spacetimes can be 
carried out. 

An alternative model is based in Chaplygin gas, whose 
equation of state is given by p = —K/p. Another possibil- 
ity is generalized Chaplygin gas (GCGV, whose equation 
of state is p = ~K/p a , < a < 1 @, 0, 1, ij. Cosmol- 
ogists became interested in this form of matter when it 
turned out to be a candidate for combining dark matter 
and dark energy: in early times Chaplygin gas behaves 
like matter and in later times like a cosmological con- 
stant. To support a wormhole, we must have [6| 



K < 



1 
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where r = ro is the throat of the wormhole, defined be- 



low. 

Two recent papers [HI EH discussed wormhole solu- 
tions that depend on a variable equation of state parame- 
ter: p/p = —K(r), where K(r) > 1 for all r. The variable 
r refers to the radial coordinate in the line clement 



ds 2 = _ e 2/(r) rft 2 



dr 2 
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dr 2 
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In this form of the line element, / = /(r) is called the red- 
shift function. The minimum radius r — ro corresponds 
to the throat of the wormhole, where b(ro) = ro. 

Ref. [12j assumes that K is both space and time de- 
pendent, i. e., K = K(r,t). Evolving equations of state 
are also discussed in Refs. [H, [13, EH- 

An interesting solution of the Tolman-Oppenheimer- 
Volkoff equation in a Chaplygin-gas setting can be found 
in Ref. [16| . Questions of stability are addressed in Ref. 

m- 

In this paper we return to the form K = K(r). The 
purpose is to show that the parameter K{r) = Br 2 ^ 1 ^ 
in the equation of state leads to a unified model for the 
generalized Chaplygin gas and phantom energy worm- 
holes, also generalized. More precisely, the equation of 
state, a motivation for which is given in the next section, 
is 



P : 



_ Br 2(a- 



and where B is a positive constant. As already noted, 
— 1 < a < corresponds to the generalized Chaplygin 
case, so that, analogously, a > generalizes the phantom- 
energy case. (Recall that the latter normally assumes 
that a — 1). 

We are striving in all cases for exact solutions. While 
such solutions are not usually obtainable, they have the 
advantage of being relatively easy to state and to an- 
alyze, without being exhaustive. On the contrary, the 
exact solutions suggest the existence of a whole class of 
solutions with a variable parameter that could model the 
properties described here. To obtain exact solutions, the 
rcdshift function 7 — 7(7") in the general line element 



ds 2 



-e 2 ^dt 2 + e 2a ^dr 2 + r 2 (d9 2 + sm 2 8dtf) (3) 
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needs to be inserted "by hand," as in Ref. [18[ | . 

From the Einstein field equations in the orthonormal 
frame, G & p — 8%T & s, the components of the Einstein 
tensor are proportional to the components of the stress- 
energy tensor. In particular, T# — p, T?f — p, and 
Tgg — Ti^ — pt, where p is the energy density, p the 
radial pressure, and pt the lateral pressure. The weak 
energy condition (WEC) is given by T & pp a p lB > for all 
time-like vectors and, by continuity, all null vectors. For 
the radial outgoing null vector (1, 1, 0, 0), the WEC now 
becomes p + p > 0. So if this condition is violated, we 
have p + p < 0. 



II. VARIABLE EQUATIONS OF STATE 

The first step in this section is to list the components of 
the Einstein tensor in the orthonormal frame [1911: 



Gfc = H e -**v (r)+ i (1 



), 
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(4) 
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G§e = Gu = e ~ Mr) (l"(r)-j'(r)a'(r) 



+h'(r)] 2 + -V(r) - -a'(r) 



(6) 



We observe next that the general line element ([3|) can 
also be put into the form of Eq. © : 



ds 2 = -e 2 ^dt 2 



dr 2 



■dr 2 



So 
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and b(r) 
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Since b(ro) — ro, ol must have a vertical asymptote at 
r = ro: lim r ^ ro+ a(r) = +oo, while 7 = 7(7") must never 
be zero to avoid an event horizon. In addition, the shape 
function b = b{r) must satisfy the flare-out conditions 
at the throat 6(r ) = r and &'(ro) < 1. Another 
requirement is asymptotic flatness: b(r)/r — > as r — > 
00. Finally, using Eqs. |(5J) and one can readily verify 
that 



b'(r) 
8nr 2 ' 



(9) 



The need for a variable K arises quite naturally: the 
first specialization, 7(7") = |ln ~, suggest ed by Zaslavskii 
[ITj], yields for a constant K > 1 dEH^I 



= 2a(r) 



1 
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an exact solution for a phantom-energy wormhole (p 
-Kp). So from Eqs. © and ©, 



-if(r) 
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Hence 



1 



b'(r) 

In particular, from Eqs. (|T0|) and 
6(r) = r 1 - (l - 1 



1 - — 



which yields b'(r) — 1/K. Since K > 1, it follows that 
b'(ro) < 1, as required. 

Using the generalized equation of state, p = —K(r)p a , 
we obtain 



K(r) 



P 



l/(87rr 2 



p a [6'(r)] a /[(87r)°r 20 ] 



{8n) a - 1 r 2a - 2 
[b'(r)] a 



So if b' = 1 /X again, then 

A'(r) = X Q (87r) Q - 1 r 2a - 2 . 



(11) 



As a consequence, we will assume from now on that the 
equation of state has the form 



P : 



where B is a positive constant. It will be seen later that 
a can be either positive or negative. 



III. THE FIRST TWO SOLUTIONS 

From Sec. |TT] the equation of state is given by 

p = -Br^ a -^p a , a ^0. (12) 
The Einstein field equations G & = 8irT & now yield 

P=^G U and ±-Grr=P=[-K(r)]p a 1 



whence 



l « 7 'W - \ (l - e- 2a(r) ) 



-1 (V**V(r) + i (l 
877 \ r r z V 



= -Br 

-2a(r) 



2(a-l) 



(13) 



Zaslavskii's function j(r) — iln £ allows us to solve 
Eq. (fT^|) by separation of variables: substituting Y(r) = 
— l/2r, we get directly 

8tt r 2 

i(^'V(.) t i(.-r>"))) 
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Raising each side to the power - produces 



1 



= B l/a r 2(l-l/o) 

r 2/a \(8n) 1 - 1 / a 
"2 



1 



x 



"Wo'W + ^ (l - e- 2Q «) 



iB-V-rsTr) 1 - 1 /" = 2e- 2a ^a'(r) + -(l- e^ r A 
r r V / 

Rearranging, we obtain 

2e- 2Q ( r V(r) 1 



and 



fl- 1 /«(8 7 r) 1 - 1 /« - (1 - e ~ 2a 



2a' {r) 



(r)) 



1 



e 2«(r) [s-i/a(87r)i-V« - 1] + 1 r ' 
The integral formula 



= u - In (Ae u + 1) 



Ae M + 1 
now yields 

2a(r) -ln{e 2a W B- 1 / a (8 7 r) 1 - 1 / a - 1 + l|= In cr 
and 

cr {e 2Q(r) [fl-V-fgTr) 1 - 1 / - l] + l} = e 2a «, 
whence 



-2a(r) _ 



1 

cr 



Since Una r ^ ro+ a(r) = +oo, e - 2a ^ = 0. Thus 



1 

ro 



B- 1 / a (87r) 1 - 1 / a - 1 



Substituting c results in 



So 



J B- 1 / a (8vr) 1 - 1 / a - 1 



,2a(r) _ 



r 



[1 - B-V a (8Try- l / a ] (1 - 



(14) 



In the special case of phantom energy, a = 1, we have 
if (r) = £?, and Eq. Q3} reduces to Eq. (fTU|) . 

NextjWe need to check the flare-out conditions at the 
throat 1]: b(ro) = ro and b'(ro) < 1. (The requirement 
of asymptotic flatness, b(r)/r — > will be dealt with 
separately.) From Eq. ([5]), 

6(r) = r {l - [l - B- 1 /°(8 7 r) 1 - 1 / a ] (l - — ) } . (15) 

Evidently, b(ro) = tq. The other requirement is 

b'(r ) = B- 1 / a {9,iz) 1 - 1 ' a < 1. (16) 



A. Generalized phantom energy 

For the generalized phantom energy case, a > 0. So rais- 
ing each side to the power —a in Eq. (TTB|) reverses the 
sense of the inequality. It follows that 



B > 



1 



(87T) 



1-a ' 



a > 0. 



(17) 



In the special case a = 1, K — B > 1, as required in the 
phantom-energy model. 

From Eqs. ([9]) and (fTB]) we can confirm that the WEC 
is violated at the throat r = ro: 



p+p = /o(l - 5r 



2(a-l) n-l-i 



1 - Br, 



o 



(a _ 1} (B-V*{*Kf-V* 



= p(i-s 1 / a (8 7 r) 1 / a - 1 ; 
i 



a-l 



< 



'('-OS^rW- 1 )" 



= 0, 



since B > l/(87r) 1 - Q , a > 0. 

While the flare-out conditions are met, the spacetime 
itself is not asymptotically flat. As a result, the wormhole 
material must be cut off at some r — A and joined to the 
exterior Schwarzschild spacetime 



2M\ , , 
1 --) dt 



dr 2 



2M 



+ r 2 (d9 2 + sin 2 6d<j) 2 ). (18) 



Such matching requires continuity of the metric. As 
noted in Ref. [2l| , since the components g^ and g^ are 
already continuous due to the spherical symmetry, one 
needs to impose continuity only on the remaining com- 
ponents at r — A: 

5tt(int) ( A ) = ftt(ext) ( A ) and 9f f (int) (A) = g?? ( cxt ) {A) 

for the interior and exterior components, respectively. 
These requirements, in turn, imply that 

Tint (A) = 7«ct(A) and b int (A) = b cxt (A). 

In particular, at r = A, 



2a(r) 



1 



1 



1 - 



b(A) 



i _ IK ' 

1 A 



So we need to determine M = hb(A), the total mass of 
the wormhole for r < A: 



M 



1 



-A 



{l- [i-S-V^Stt) 1 - 1 /"] (l-^)} 



Since 7(7") = |ln -, we have 



ln (c/r) 



so that 



.4 



2Af 
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After solving for c, we find that 
f 



- {a- a (i - [1 - b- 1 ^) 1 - 1 / ] (i - 5)) } 

for r < r < A. At r = A, e 2 ^ _ ! _ /A, as 
required. So for r > A, the metric becomes 



ds 2 



6(A) 



6(A) 



dr 2 



+ r 2 (d6» 2 + sin 2 6ld0 2 ). 



B. Generalized Chaplygin gas 

For the Chaplygin gas case we have 

p= -Br 2{a - 1) p a , a<0. 

To determine the condition on B, it is best to return to 
Eq. UJJ and observe that 



B- 1 '" 1 ^) 1 - 1 ^ < 1. 



(19) 



Since —a > 0, raising each side to the power —a retains 
the sense of the inequality. The result is 



B < ^ , a < 0. 

(8^^) 1 - a, 



(20) 



This condition is similar to K < (87T?' 2 ,) ( a+1 ), a > 0, 
in the original GCG model. In fact, if B is replaced by 
B*r Q 2( * a l ' in the equation of state, then from Eq. (fl^j) . 



B* < 



„2\l-a ■ 



a < 0. 



It also follows from inequality (f2TJ)) that the WEC is 
violated at the throat (Sec. IIII A[) . 

Conversely, ifl/K = B- 1 / a (8n) 1 - 1 / a in Eq. QU, then 
B = K a (8n) a - 1 , which agrees with Eq. (jTT]) . 



IV. THE NEXT TWO SOLUTIONS 

As noted earlier, there are only two ways to insert j(r) 
by hand. The second is to assume that 7' = 0. If K is 
constant, this yields Lobo's solution \^ 



1 



1-1/A' ' 



(21) 



Returning to Eq. (|13p . we now have 
i_ (1 _ e -a«(r)) ^Sr 2 ( a - 1 )(8 7 r) 1 - a 

'e- 2a < r V(r) + 1 (l - e~ 2Q M) 



(22) 



Proceding as before, we can put Eq. ([22)) into the fol- 
lowing form: 



2e- 2a ^a'(r) 

I _ e -2a(r) 



-l + B- 1 /a(8 7r )i-i/«(l- e -2«M) 1 / a 1 r' 

Integrating the left side requires the formula 
du 1 



(23) 



u{-l + Au B ) B 
The result is 



[ln(Au B - 1) - In 71 s ] 



B 



By letting r — r , we obtain 
1 



ro 



- 1 



a/(l-a) 



and, finally, 



where 



2a(r) _ 



1 - [F(r)} 



a/(a-l) ' 



(24) 



F(r) = 



(?) 



(a-l)/a 



1) 



(a-l)/a 



In the generalized phantom energy case observe that solu- 
tion (|24|) looks substantially different from solution (|2ip . 
It can be shown, however, that 



lim e 2a{r) = - ■ 



(25) 



(recall that B = K when a = 1.) The reason for this 
outcome can be seen by means of the following heuris- 
tic argument: in Eq. (|23|) . for a very close to unity, let 
H = (1 — e -2a(r)^ 1 /i-i ^ go that q~ 1 Treating fiasa 
constant, Eq. (|2"3")l can be integrated to yield 



2a(r) 



1 



l-S-V«(87r) 1 - 1 /afJ ' 



1 



(26) 



If a = 1, so that SI = 1, we obtain solution (|2"5 
From 



6(r) =r Ib- 1 ' 11 ^) 1 - 1 ^ 



ro^(a-l)/« 

r 



(7) (27) 
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one can show that 

&(r„) = r and 6'(r ) = B^' 11 ^) 1 - 1 ' 11 , 
as before [Eq. (fT6|)]. So b'(r Q ) < 1 provided that 

B > — — !- — , a > 0. 

(Stt) 1 - 

For the corresponding Chaplygin case, 

B < 7— 4i — , a < 0. 

(Stt) 1 - 

As a consequence, the WEC is violated at the throat in 
both cases, already shown in Sec. IIIIl 

For the junction at r = A, we determine b(A) = 2M, 
as before. Then the metric becomes 



dt 2 + e 2a ^dr 2 



+ r 2 (dd 2 + sin 2 8 dej) 2 ) 



for ro < r < A, and 



ds 2 



6(A) 



r 

r 2 (d8 2 + sin 2 8 d(j) 2 ) 



dr 2 

2/ 



for r > A. 

Conversely, suppose b(r) has the form 



CONCLUSION 



This paper introduces a new equation of state with a 
variable parameter: 



p= -Br 2{a - 1] p a , a^O, B>0. 
This equation of state generalizes the notion of phan- 
tom energy. The result is a combination of the general- 
ized phantom energy and the generalized Chaplygin gas 
models (themselves generalized) to describe a traversable 
wormhole geometry. The traversability conditions can be 
determined from the WEC violation in Sec. IHI Al 



p fi _ B^ a {^f' a - l \ < 



whenever B 1 ^ > l/(87r) 1 / a 1 . This inequality implies 
the following: if 



B > 



then a > 0, 



and if 



b{r) = r 



1 

K 



1 - 



(a-l)/a 



(t) 



based on Eq. (QTJ). From Eq. © 
1 f 1 



V - 



8-Kr 2 1 K 



i iw n a/ (a— 1) 



Then, using p = b' (r) / (8irr 2 ) , the expression for 
—K{r) = p/p a can be reduced to the form 

K(r) ^ K a (8?:) a - 1 r 2a - 2 , 

in agreement with Eq. (jlip . 



in) 1 - 11 



then a < 0. 



In the first case the wormhole is supported by general- 
ized phantom energy and in the second by generalized 
Chaplygin gas. 

The solutions obtained are exact, yielding easily stated 
conclusions, but extracted at a price: the redshift func- 
tion must be inserted "by hand." Yet apart from the 
ability to find exact solutions, there does not appear to 
be anything artificial or even special about these redshift 
functions. Instead of producing purely mathematical so- 
lutions, it seems more likely that the generalized phan- 
tom energy and Chaplygin gas models can, at least in 
principle, be unified by equations of state having a vari- 
able parameter. 



[1] M.S. Morris and K.S. Thorne, Am. J. Phys. 56, 395 [3] S.J. Perlmutter et al, Astrophys. J. 517, 565 (1999). 

(1988). [4] R. Bousso, arXiv: 0708.4231. 

[2] A.G. Riess et al, Astron. J. 116, 1009 (1998). [5] S.V. Sushkov, Phys. Rev. D 71, 043520 (2005). 



6 



[6] F.S.N. Lobo, Phys. Rev. D 73, 064028 (2006). 

[7] A.Y. Kamenshchik, U. Moschella, and V. Pasquier, Phys. 

Lett. B 511, 265 (2001). 
[8] M.C. Bento, O. Bertolami, and A. A. Sen, Phys. Rev. D 

66, 043507 (2002). 
[9] P.K.F. Kuhfittig, arXiv: 0802.3656. 
[10] F.S.N. Lobo, Class. Quant um Gray. 24, 2 401 (2007). 
[11] F. Rahaman et ai, arXiv: |gr-qc/0701032j 
[12] P.K.F. Kuhfittig, Schol. Res. Exch. 2008, 296158, (2008). 
[13] U. Mukhopadhyay, S. Ray, and S.B. Dutta Choudhury, 

arXiv: 0711.4800. 
[14] J. Sola and H. Stefancic, J. Phys. A 39, 6753 (2006). 
[15] J. Sola and H. Stefancic, Mod. Phys. Lett. A 21, 479 



(2006) . 

[16] V. Gorini, U. Moschella, A.Y. Kamenshchik, V. Pasquier, 
and A. A. Starobinsky, Phys. Rev. D 78, 064064 (2008). 
[17] E.F. Eiroa and C. Simeone, Phys. Rev. D 76, 024021 

(2007) . 

[18] P.K.F. Kuhfittig, Class. Quantum Grav. 23, 5853 (2006). 
[19] P.K.F. Kuhfittig, Phys. Rev. D 66, 024015 (2002). 
[20] O.B. Zaslavskii, Phys. Rev. D 72, 061303(R) (2005). 
[21] J.P.S. Lemos, F.S.N. Lobo, and S.Q. de Oliveira, Phys. 

Rev. D 68, 064004 (2003). 
[22] F.S.N. Lobo, Phys. Rev. D 71, 084011 (2005). 



